Measurement of two-qubit states by quantum point contacts 
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We solve the master equations of two charged qubits measured by two serially coupled quantum 
point contacts (QPCs). We describe two-qubit dynamics by comparing entangled states with product 
states, and show that the QPC current can be used for reading out results of quantum calculations 
and providing evidences of two-qubit entanglement. We also calculate the concurrence of the two 
qubits as a function of dephasing rate that originates from the measurement. We conclude that 
coupled charge qubits can be effectively detected by a QPC-based detector. 



I. INTRODUCTION 

Quantum information processing in charge-based solid 
state nanostructures has attracted wide spread attention 
because of the potential scalability of such devices, and 
the relative ease with which such charge devices can be 
manipulated and detectedi 1 ' 2 ' 3 ' 4 Recently, two-qubit co- 
herent evolution and possibly entanglement have been 
observed in capacitively coupled Cooper pair boxes A For 
universal quantum computing, two-qubit operations are 
required, so that the realization of controlled two-qubit 
entanglement is regarded as a crucial milestone for the 
study of solid state quantum computing. While two- 
qubit information can be detected with one measurement 
device on each qubit, it is also important to search for a 
detector that is directly sensitive to two-qubit informa- 
tion, and to develop a proper formalism to study two- 
qubit measurement processes AL^ 

The ultimate criterion for the detection of qubits is 
whether we can distinguish the results of a quantum com- 
putation by the output signal of the detector, e.g. cur- 
rent or conductance of a single electron transistor. In 
the case of one qubit, two single-qubit states | |) and 
| |) need to be distinguished. In the case of two qubits, 
four two-qubit state, | ||), | It): I 11), an( i I TT) ( we wm 
call them \A) ~ \D}) need to be distinguished before the 
qubit states are destroyed by the measurement. As we 
mentioned above, measurement of multi-qubit states can 
generally be achieved by multiple single-qubit measure- 
ments on each of the qubits, respectively. Here we study 
a different detection process: the temporal behavior of 
a detector (QPC in the present study) that simultane- 
ously couples to two qubits. We show that information 
contained in the temporal evolution of the QPC current 
can help us distinguish different two-qubit product states, 
and some entangled states from the product states. In- 
deed, one motivation of our study is to obtain direct ev- 
idence for the entanglement of the qubits, possibly from 
the detector current or other measurable quantities. 

In Ref. 0, we studied a particular scheme for the 
quantum measurement of two charge qubits by a two- 
island single-electron transistor (SET), and showed that 
the SET is an effective detector of the two-charge-qubit 
states. Here the charge qubits are constituted of two cou- 
pled quantum dots (QD) with one excess electron. Due 



FIG. 1: Two charge qubits (using double quantum dot 
charged states) are capacitively coupled to a detector of two 
serially coupled QPCs. J is the strength of inter-qubit inter- 
action. No tunneling is allowed between the QPC detector 
and any of the qubits. 



to tunnel coupling of the QDs, the wave functions in a 
qubit are superpositions of localized states from each of 
the QDs. If a qubit is prepared in a single QD state, 
it tends to oscillate between the two sides of the dou- 
ble QD. If we define the local states as | f) and | J.), 
the qubit state oscillates between the two logical states 
with a frequency determined by the tunneling coupling 
and the difference of the energy-levels of the two QDs. 
Time-dependent behavior of this coherent oscillation of 
the qubits is determined by the initial state. If we take 
the initial time to be that when a final quantum opera- 
tion is applied to the qubit, the detector readout current 
reflects the results of quantum calculation. The qubit 
states interact with the readout current by changing the 
energy (and therefore occupation) of the electronic states 
in the SET islands and possibly the tunneling rates of the 
junctions (by modifying the island electronic states them- 
selves) on both sides of the islands. Although, in Ref. 0, 
we show that the SET can distinguish the different co- 
herent oscillations between the two-qubit product states 
and the entangled states, we have not yet investigated 
the two-qubit dynamics itself. 

Here we would like to study quantum point con- 
tact (QPC) as a detector for two coupled charge qubits 
(Fig.QJ. A QPC is essentially a one-dimensional conduct- 
ing channel and is considered to be an effective charge de- 
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FIG. 2: Electronic states in the inter-QPC island. We assume 
that there is only one accessible electronic orbital state on the 
island. There are thus totally four possible island states: (a) 
Empty dot — state "a" has no excess electron on the dot. (b) 
Single-electron dot — state "6" has one electron and is spin- 
degenerate, (c) Two-electron dot — state " c" has two electrons 
in a spin singlet state occupying the same orbital state. 



tector, similar to the SET. The particular scheme we con- 
sider consists of two low transparency QPCs connected 
in series through a single level quantum dot. Each of 
the QPCs is close to a charge qubit so that its cur- 
rent/conductance is dependent on the state of the re- 
spective qubit. Compared to the SET detector, the QPC 
detector interacts with qubits only through the change 
of tunneling rates. Although the SET detector is able 
to describe a variety of features of the internal states of 
qubits in Ref£, we could not identify which of the two 
interactions (that between the islands and qubits, which 
modify the level occupations on the SET islands; or that 
between the tunneling junctions and qubits, where po- 
tential by the qubits modifies the tunneling rates) plays 
the major role in the SET detector. Thus, an impor- 
tant question is whether or not the QPC detector that 
interacts with qubits only by the change of the tunneling 
rates is also an effective apparatus for detecting the qubit 
states. In this paper, after discussing the basic two-qubit 
dynamics with no detector, we focus mainly on the fol- 
lowing issues: (1) whether we can distinguish the four 
product states \A) ~ \D) of two coupled charge qubits in 
the time-domain with a serially coupled QPC detector, 
(2) whether we can distinguish the entangled states from 
the product states of these two qubits, and (3) whether 
the quantum Zeno effects in the coupled charge qubits 
can be observed. In the following sections, we solve the 
master equations for the coupled qubit-QPC system and 
investigate the effectiveness of the proposed QPC detec- 
tor. In Sec. II, we present our formulation of two qubits 
and the QPC detector. In Sec. Ill, we discuss the dif- 
ference between a dynamics of single qubit and that of 
two qubits. In Sec. IV, we show the numerical results of 
two-qubit detection by QPC. Sec. V is devoted to dis- 
cussion for the QPC detection, and Sec. VI consists of a 
conclusion. 



II. FORMULATION 

In the present measurement scheme, the QPCs are 
capacitively coupled to the charge qubits (Fig. 0), so 
that the current through them sensitively depends on 
the states of the qubits. We describe the two QPCs us- 
ing two tunnel matrix elements only and neglect further 



internal structures AiS We assume that the qubit-QPC 
coupling is purely capacitive, so that there is no current 
flowing from the qubits to either of the QPC electrodes. 
The Hamiltonian for the combined two qubits and the 
QPCs is written as H = £f q b+-ffqp C +-ffint- -Hqb describes 
the two interacting qubits (left and right, as illustrated 
in Fig. [I}, each consisting of two tunnel-coupled QDs and 
one excess charged 

) + JvLz<JRz, (1) 

a=L,R 

where Q a and A a are the inter-QD tunnel coupling and 
energy difference (gate bias) within each qubit. Here we 
use the spin notation such that a ax = a\jj a + b\a a and 
a az = a\a a ~ b\b a (a = L, R), where a a and b a are the 
annihilation operators of an electron in the upper and 
lower QDs of each qubit. Thus, | f) and | J.) refer to 
the two single-qubit states in which the excess charge is 
localized in the upper and lower dot, respectively. J is 
a coupling constant between the two qubits, originating 
from capacitive couplings in the QD system^ 

The two serially coupled QPCs are described by H qpc : 

ct = L ,R 1 

+ E d4d s + Ud\d ] d\d l ■ ( 2 ) 

s=T4 

Here Ci LS (ci RS ) is the annihilation operator of an elec- 
tron in the i^th (i/jth) level (iz,(*ii) = 1, n) of the left 
(right) electrode, d s is the electron annihilation operator 
of the island between the QPCs, Ei LS (Ei RS ) is the energy 
level of electrons in the left (right) electrode, and Ed is 
that of the island. Here we assume only one energy level 
on the island between the two QPCs, with spin degen- 
eracy. Vi LS {Vi RS ) is the tunneling strength of electrons 
between the left (right) electrode state ils (irs) and the 
island state. U is the on-site Coulomb energy of dou- 
ble occupancy in the island. Note that the number of 
island states here (Fig. |2J) is much smaller than that of 
the two-island SET states^ where we need to take into 
account at least 10 types of detector states. In Ref. 0, we 
observed that the two-island SET can represent a variety 
of qubit states because of its large number of degrees of 
freedom. With a much simpler state structure for the 
present coupled QPC scheme, we will study whether the 
QPC current could still represent the qubit states with 
sufficient clarity. 

The capacitive interaction between the qubits and the 
QPCs is represented by H- mtl which contains the infor- 
mation that localized charge near the QPCs increases 
the energy of the system electrostatically, thus affecting 
the tunnel coupling between the QPCs and the island in 
between^ 

H '^= Y, Y. 8V ^ c Ls d s + dlc laS )a az , (3) 

a—L,R i a ,s 
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where SVi aS is an effective change of the tunneling 
strength between the electrodes and QPC island (we 
shift the origin of the interaction energy such that 
SVi a , s — corresponds to the case where qubits are 
not polarized a az — 0). We assume that the tun- 
neling strength of electrons weakly depends on the en- 
ergy Vi aS = V a (Ei aS ) and electrodes are degenerate up 
to the Fermi surface. Then the qubit-QPC interaction 
dictates that qubit states influence the QPC tunneling 
rate in the form of T^\e) = 2np a (E)\ F Q (±) {E)\ 2 and 

Ti ±Y (E) = 2^p a {E+U)\vt\E+U)\ 2 , where V^\e) = 
V a {E) ±6V a (E) {SV a (E) > 0), and p a (E) is the density 
of states of the electrodes (a — L,R). Hereafter, we use 

la 's and s estimated at the Fermi energy /i Q of the 
electrodes to describe the tunneling rate in the detection 
process of the qubit states by the two QPCs. This is 
reasonable from a practical standpoint since many ex- 
periments are described using Top*^- The values of the 
corresponding ri^s are determined by the structure of 
the system such as the distance between the qubits and 
the QPCs. For example, a 1 1) state (| |) state) in a qubit 
means the excess charge is localized in the lower (upper) 
dot, so that the corresponding tunneling rate should be 
r L } (^l )- Therefore, two-qubit state \A) would lead to 
QPC tunneling rates of (r^\ r^ } ) or (r^' , r^'V de- 
pending on whether or not the island between the QPCs 
is doubly occupied. Similarly, \B), \C) and \ D) states cor- 

respond to (r«, if) [or (p"', if')], (if, if) 
[or(lf Py')],and(r«, if) [or (if if')], re- 



spectively. 

We construct the equations for the qubit-QPC density 
matrix elements at zero temperature T — 0, following the 
procedure developed by Gurvitzi2iii This method is ap- 
plicable when the energy level of the inter-QPC island is 
in between the chemical potentials of the two electrodes. 
The wave function IvE'(t)) of the qubit-QPC system can 
be expanded over the Hilbcrt space spanned by the two- 
electron states of the qubits, the island states of the QPC 
shown in Fig. ^ an d all possible electrode states. We 
choose |0) to refer to the initial ground state of the whole 
detector system (two electrodes and the inter-QPC is- 
land) where the two electrodes are filled with electrons 
up to /ii and jitR, respectively, and the two QPCs and 
the inter-QPC island are empty of excess electrons. The 
basis states for the QPC can then be constructed from |0) 
by moving electrons from the left electrode (with higher 
chemical potential) to the inter-QPC island and the right 
electrode. We categorize the detector states by the num- 
ber of electrons that are transferred from the left to the 
right electrode (Fig. 1): 

|*(t))=|*o(t)> + |*i(t)>, (4) 

where [^(i)) is the part of the wave function that 
no electron tunnels through to the right electrode and 
|\I/i(i)) represents the part of the wave function where 
one or more electrons are transferred to the right elec- 
trode. |^o(t)) can be expressed as 



i*o(<)) = E { 6 (0)a ' 2 w + E & L 0)fc ' 2 w d ^+E 6 LtuW 44^4 io>i*>, (5) 

z=A,B,C,D I Is hh ) 

where M°) a,z (i), (t) and bf^J^(t) are coefficients for the respective states. The superscripts refer to the number 

of electrons transferred (0 here), the states of the QPC island (as illustrated in Fig. |2J), and the four two-qubit basis 
states. The subscripts refer to the left electrode states from which electrons tunnel into the islands. Thus each of the 
terms in |\Po(i)) indicates a state with as little as zero but up to 2 electrons moved from the left electrode to the QPC 
island, while no electron is transferred to the right electrode. can be expressed as 

oo ( \ n 

i*i(*)> = E E Siw+E&w ^ s +EfeV. /3 „w d\d\c hr c l2l ®n (4 s , Cr v) io>]^> , (6) 

n=l z=A,..,D y is i x i 2 ) i = i 

01 --Pn 



where = (Z,-,r^,s-) represent the left electrode, right 
electrode, and spin states involved in the transferred elec- 
trons. Similar to the expressions of the coefficients for 
|* (t)>, here b%>X(t), b££'(t) and b^^t) are 
coefficients for the states with n electrons transferred to 
the right electrode, and another to 2 electrons moved 
from the left electrode to the QPC island. The super- 



scripts again refer to the number of transferred electrons 
(n), the QPC island states, and the qubit basis states. 

Substituting this wave function into the Schrodinger 
equation for the whole qubit-QPC system, i\^(t)) = 
H\ty(t)) (having set h = 1), we obtain a set of alge- 
braic equations for the coefficients in Eq. (J5J and Eq. 
Assuming zero magnetic field and spin-independent 
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electron tunneling, the density matrix elements can be 
defined as 

J 

(7) 

where b^ 3 1 ' zl {E) is a Laplace-transformed element of 
b W p' zl {t) after summing over (3 — {0, /3i, /?2, • • • , ■ • •}, 
the electrode states of transferred electrons as discussed 
above ('0' corresponds to coefficients in Eq.JSJ)). 

After a lengthy calculation, we obtain 48 equations 
for the density matrix elements p^ 1Z2 (t), where u — 
a,bi,bi,c indicate quantum states of the inter-QPC is- 
land (Fig. I2J) as shown in Appendix A. 12 . Because we as- 
sume that there is no magnetic field, p z \ Z2 {t) — p Zl * 2 {t). 

The readout current is obtained as a time derivative of 
the number of electrons in the island&tLil 

J(t) = eT R [p b] {t) + p H {t)} + 2eT' R p c {t) , (8) 

where p u given by p u = p AA + Pu B + Pu C + Pu° is thc 
occupation probability of the particular island state it. 

As the difference between and Ta ' increases, the 
current difference that depends on the difference of qubit 
states increases as well, while the qubits lose their coher- 
ence faster due to the fluctuations in the QPC current, 
which lead to fluctuations in the qubit energy levels and 
thus dephasing. We quantify the strength of the mea- 
surement by dephasing rates defined as 

If ee (v^-vtF) , (9) 

where a = L,R. These rates are the coefficients of 
the off-diagonal density-matrix elements of the time- 
dependent reduced density matrix equations for the 
qubits. The reduced density matrix elements are 

p ziz2 =pr 2 +pir+ptr+pc Z2 - ( 10 ) 

This definition of dephasing rate is originally introduced 
by Gurvita^ for the case of a single qubit. The dephasing 
time is taken as coinciding with the measurement time. 
Compared with Ref.|9|, where there is a single off-diagonal 
density- matrix element, we cannot define a single dephas- 
ing rate because of the complexity of our density-matrix 
equations. 

Current of a QPC in the tunneling regime is very sen- 
sitive to the potential variations of the QDs that are set 
close to the QPC channel^ We thus can choose from a 
wide range of parameters for our QPCs. Here we use a 
particular set of representative parameters: 

= r£ = r2=rg=r(-)=o.sr, (11) 
r c = r£ = rg=rg=r< + )=i.2r, (12) 



which lead to = Tf(= T d ) ~ 0.04r for a typical 
case (r is a tunneling rate in the absence of the qubits, 
so that dephasing rate is more than one order of mag- 
nitude smaller, corresponding to a relatively weak mea- 
surement). We can regard T^" 1 as the typical measure- 
ment time. Obviously, the qubit dynamics (represented 
by tunneling rate £1) would be able to be detected when 
Q^ 1 < T^ 1 . Because, in the present setup, the cur- 
rent without qubits saturates in the range of ~ 
the time T^ 1 would serve as a standard of when mea- 
surement started. We can also incorporate the effect 
of Coulomb interaction by setting = as a lim- 

iting case of strong on-site Coulomb blockade (U — ► oo 
in Eq. (0 so that no double occupation is possible) , while 
for weak Coulomb interaction on the island we can set 
r£ ±) = ri ±] at the limit of U = 0. 



III. QUBIT DYNAMICS WITHOUT DETECTOR 

In order to better understand our numerical results 
and the capability of our QPC detector, it is instructive 
to first examine the dynamics of both a single qubit and 
two qubits in the absence of any detector, and discuss 
how the qubit dynamics is measured by the detector. 

We first solve the density matrix equations for a single 
qubit on the basis of localized single quantum dot states 
I T) and | I): 

Pit = «fi(pn -AlT*K (13) 
PU = ^(PlT -Pu)> ( 14 ) 
PU = iA PU +^(PTT "Pu)- ( 15 ) 

For the simple case of A = (no voltage bias between the 
two dots so that qubit dynamics is completely determined 
by the inter-dot tunnel coupling Q, which corresponds 
to the optimal operational point in terms of minimum 
dephasing as discussed in Ref«ii), and starting from one 
of the localized states "["-state (ptt(t = 0) = 1) or J.-state 
(Pu(* = 0) = 1), we have 

p TT (i) = p TT (0)cos 2 (ffi)+p u (0)sin 2 (m), (16) 
PllW = Pu(0)cos 2 (m)+p TT (0)sin 2 (f!i), (17) 

PU(*) = J(PTt(°) -PU (0))sin(2fit). (18) 

These solutions dictate that the QPC current should es- 
sentially be determined by p-\i(t) — Pu(t) = [pfr(0) — 
pjj (0)] cos 2f2£. The oscillatory component of the QPC 
current should thus be d ominated b y a 2fl component 
(in the case of A ^ 0, 2y / fl 2 + A 2 /4), and the temporal 
evolution of the current is intimately related to the initial 
state. 

We can also infer information on the two-qubit prod- 
uct states from the detector current in a similar manner 
because density matrices of the product states are writ- 
ten as p AA (t) — pf, (t)pf\(t) and so on. Here we solve the 
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two-qubit dynamics in the absence of the detector by ex- 
panding the density matrix on the basis spanned by the 
Bell states: |ei)=(|U) + | Tt»/V2, |e 2 )=(|U> - 1 TT»/V2, 
|e 3 ) = (|IT) + in))/V2, and |e 4 ) = (| 4?) - I U»/V2 (sin- 
glet state). If we assume two identical qubits (J7l = £Ir 
and Al — A^(— A)), the density matrix equations for 
the two qubits (without the QPC detector: 1^ = 0) are 
written as 

p eiei = 

p e2fi2 = 2zA(p e2£l - p e3£2 ) (19) 

p e 2 e 4 = _ 2i J p e2C4 _ 2iA P ^ ei 

p eiei = 2itt(p eie3 - p e36l ) + 2iA(p eie2 - p e2Cl ) 
pw* = -2iQ(p e ^ - p e3fil ) (20) 

p e ie3 = _2iO(p e363 - p eiei )-2iJp eie3 -2iAp e2 e 3 



p eie 2 


= -2iSlp e3B2 


-2iA(p eiei -p e2e2 ) 




= 2iflp e ^ - 


2iJp e ^-2iAp e ^ 




= -2iVLp eiei 




p eie-i 


= -2iflp e3e4 


-2iJp eiei ~2iAp eie - 1 



If A = 0, which again corresponds to the optimal oper- 
ational point, the density matrix equations can be clas- 
sified into four groups (indicated by the four parenthe- 
ses above). First of all, Eqs. I|19fl shows that the prob- 
abilities in \e2) and \ei) states are time-independent. 
On the other hand, according to Eq. H20J1 . the proba- 
bilities in \ei) and |es) states oscillate as a function of 
{cos(4fi*i),sin(4fi*t)} (fi* = y/Q 2 + J 2 /A). Meanwhile, 
Eqs. I|21|l and IL'l't indicate that the off-diagonal elements 
{p eie2 ,p e2E3 ,p e3e4 ,p eie4 } contain {cos(2fi*t),sin(20*i)} 
type of oscillations. Therefore, the occupation probabil- 
ities for the product states, p AA = (p eiei + p eiC2 + p eie2 + 
p e2£2 )/2, p BB = (/Z 363 +p e3Ci +p e3Ci + p eiei )/2, p cc = 
(p e3e3 -p e3e4 -p e3e4 +/9 e4e4 )/2, and p DD =(p e ^ - p^ 2 - 
p eie2 +p e2£2 )/2, should all contain {cos(2ft*i), sin(2fi*i)} 
oscillations, reconfirming the calculations on single-qubit 
dynamics at the beginning of this section. Therefore, 
we should be able to distinguish pure entangled states 
from pure product states \A) ~ \D) based on whether 
the detected period of the coherent oscillations is lim- 
ited to {cos(4fi*i), sin(4fi*t)} (|ei) and |e 3 )) or time- 
independent (\e2) and \ei)) in the limit of weak inter- 
action between the qubits and the QPCs. Such behavior 
is indeed evident in our results as shown in the following 
section. 



IV. NUMERICAL RESULTS OF QPC 
DETECTION 

In Ref0, we clarified three major issues regarding the 
capability of the two-island SET by monitoring its time- 
dependent readout current: (1) the two-qubit product 
states (eigenstates in the absence of inter-qubit interac- 
tion and inter-dot coupling within each qubit) \A) ~ \D) 
can be distinguished; (2) pure entangled states and pure 
product states can be distinguished; (3) quantum Zeno 



effect is present in a two-qubit system. In the following 
we show that similar results are obtained for the serially- 
coupled QPC detector despite its simpler state structure. 

Figure [3] shows the time-dependent current at small 
time t ~ in the case of weak Coulomb interaction 
(U = 0) (ri^ = Ta) assuming that initially the two 
qubits are in one of the four product states. To calcu- 
late the current when the two-qubit initial state is \A), 
for example, we set &( ) a '^(0) = 1 and the other coeffi- 
cients to zero in the total wave function (Eqs. (JjJJ and 
©), which means that p AA (0) = 1 and other density 
matrix elements are all zero at t = 0. At small t ini- 
tial state \A) (with both electrons located in the respec- 
tive lower dots) leads to the strongest suppression of the 
QPC current, while initial state \D) (with both electrons 
located in the respective upper dots) the least. States 
\B) and \C) also produce different QPC currents. The 
reason is that there is a finite bias between the left and 
right electrodes, so that current flows only in one direc- 
tion. Consequently, |C), with the left qubit electron in 
the upper dot (thus less suppression on current), pro- 
duces a faster rise in current than \B). Since the product 
states are not the two-qubit eigenstates, they evolve into 
superposition states and the corresponding QPC current 
oscillates. Nevertheless, we can distinguish the four ini- 
tial product states by the values of the readout current 
in both J = and J ^ cases. Hereafter we will fo- 
cus on the J = case. As shown in Fig. [21 the current 
differences between the four two-qubit states can be de- 
tected before the coherent motion of the qubits changes 
the two-qubit state as fit < ir /4. 

One observation we made for charge qubits measured 
by an SET detector is that the amplitude of the SET 
current oscillations corresponding to the pure entangled 
states are smaller than those of the pure product states £ 
Similar effects are also observed for the QPC detectors 
here, as indicated in Figs. l4l and II VI A qualitative reason 
is that the wave functions of the entangled states in the 
charge qubits extend over both qubits compared to the 
product states, so that the charge distribution of the en- 
tangled states is less effective in influencing the readout 
current. Quantitatively, for instance, Eq. (|19fl also dic- 
tates that current corresponding to a singlet state should 
have very weak time dependence. Indeed, Fig. 0] shows 
strong differences between QPC currents for the singlet 
state \e4/ and product state \B): the detector current on 
the singlet state is essentially a smoothly rising function 
of time, while the current for the product state has an 
oscillatory component of frequency ~ 2f2 at V g = 0. We 
obtained similar current behaviors for other entangled 
states and product states, where the peaks of the coher- 
ent oscillations in the other product states are shifted as 
inferred from Fig. [31 These features hold as long as the 
entangled states are close to any of the Bell states, |ei), 
~ | e^). Figure llVl shows the time-dependent current of 
the generalized singlet state cos 6\ l*\)+e lip sin 6\ tj) in the 
range of <p — n, = 6 = n/2. We found that the unifor- 
mity of the readout current holds approximately up to 
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FIG. 3: Time-dependent QPC current I(t) of the (7 = case 
(r„ =ri ± - ) ) starting from four product qubit states \A) ~ 
\D) at time t = 0. fi L = = 0.75r, r d = 0.04r. The 
two panels refer to two different inter-qubit interaction: (a) 
J = 0, (b) J = T. We can distinguish the four product states 
in both the J = case and the J = P case. This shows 
that we can distinguish the four two-qubit product states in 
a range of inter-qubit coupling strength. 




FIG. 4: Time evolution of QPC current I(i) corresponding to 
the product \B) state (panel (a)) and the entangled singlet 
state \e&) (panel (b)) when the qubit gate-bias V g (= Al = 
Ar) changes. The relevant parameters are chosen as Ql — 
Q.r = 0.75V, J = 0, U = and F d = 0.04F. The I(t) for the 
product state ((a)) explicitly reflects the coherent oscillations 
of the qubit states, while those for the entangled state are 
rather uniform. 



|0±7r/4| =tt/12, which is remarkably robust (similar to 
the case of charge qubits measured by SET detectors^). 
In addition, in Fig. IIVI the current for \C) state, another 
product state, also contains an oscillatory component of 
frequency ~ 2f2. 

An interesting aspect in studying quantum measure- 
ment is to explore the backaction of the measurement 
apparatus on the qubits. In this context, the quantum 
Zeno effect refers to the phenomenon that a continuous 
measurement slows down transitions between qubit quan- 
tum states due to the collapse of the qubit wave function 
into observed states. This phenomenon might be useful 
in quantum computation because it preserves the results 
of quantum calculations for a longer period of time. 15 
Figure demonstrates the quantum Zeno effect for two 
qubits measured by the QPC, where the initial state is 



chosen to be \D) state (p DD (t = 0) = 1). As the measure- 
ment strength increases (Yd increases), the oscillations of 
the density matrix elements of the two qubits are delayed, 
which is clear evidence of the slowdown described by the 
Zeno effect. 

In general, increasing measurement strength (i.e. the 
coupling strength between the qubits and the QPCs) 
leads to faster entanglement between the qubits and the 
measuring apparatus, so that measurement leads to pro- 
jection of qubit states into product states. Therefore, 
stronger measurement strength destroys entangled qubit 
states more rapidly. This is in contrast to the product 
states, for which the quantum Zeno effect is observed 
(Fig. EJlii 5 - We use the concept of concurrence^ to quan- 
tify two-qubit entanglement and calculate concurrence 
in the presence of increasing measurement strength. Fig- 
ure IIVI shows the effect of measurement on the singlet 
state, demonstrating that stronger measurement (in the 
form of larger Yd) degrades the entanglement (in terms 
of concurrence) more rapidly. As seen from Eas. (|13|l - 115|) 
and from (|I9|) - (|22|I . product states and entangled states 
discussed here are generally not two-qubit eigenstates 
even in the absence of the detector, and thus could evolve 
into each other through the time-dependent coherent os- 
cillations. Strong detection enhances the mixing of these 
states and makes it more difficult to infer the qubit states 
from the detector current. Figures|Hl(a) and (b) show the 
time-dependent currents of \ei) (singlet state) and |es) 
state as functions of increasing measurement strength. 
Without the detector, singlet \e4) state should be time- 
independent according to Ea. H19|) . and {es) should show 
4f2 oscillation according to Eq. (|20|l . Figures |H| (a) and 
(b) indicate these characteristics in the weak measure- 
ment case Yd < 0.04r, which is also the case that we 
discussed concerning Fig. 01 I n this region, we would 
be able to distinguish the different behaviors of entan- 
gled states and product states. However, as the strength 
of measurement increases, the detector current starts to 
acquire other oscillatory components, which means that 
both states are mixing with other states after t = 0. Fig- 
ure |S] (c) is a time-dependent diagonal matrix element 
^e 4 e 4 Q f -(.j^g sm gi e t state. This figure also shows that the 
singlet state mixes with other states as the strength of 
measurement increases. 

In the case of a strong Coulomb interaction so that 
= 0, we have obtained similar results, except that 
the magnitude of the average current is reduced by half 
because the onsite Coulomb interaction closes one trans- 
mission channel. This is different from the coupled SET 
detector we studied before, where current uniformity in 
finite-C/ model is more persistent than in the infinite-C/ 
model, because the internal degrees of freedom in the 
two-island SET allow a redistribution of electrons be- 
tween the islands. Here there is only one island with three 
island states (unoccupied, singly occupied, and doubly 
occupied, shown in Fig. EJ) . The much simpler internal 
dynamics of these states is insufficient to cause any large 
change in the QPC current when Coulomb interaction is 
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FIG. 5: A contour plot of the time evolution of QPC current 
for states ranging between \B) state and |C) state through 
singlet state |e4) (see text). The current for the "general" 
singlet state shows uniform characteristics when it is close to 
the exact singlet state |e4) as \0±tt/4\ = tv/12. The chosen 
parameters are similar to what we have before: Q.l = = 
0.75r, J = 0, V g = 0, U = 0, and r d = 0.04r. In addition, 
the current for |C) has a oscillatory component of frequency 
20 = 1.5. 



accounted for. 



V. DISCUSSION 

In our study so far we have demonstrated that two- 
charge-qubit state information can be clearly revealed 
by the transient current variations in a serially coupled 
QPC charge detector. An important question is then 
whether such current evolutions are experimentally ob- 
servable. In our calculation, T, the QPC tunneling rate, 
is the physical quantity that can be directly connected to 
experiments. For example, for T in the order of 100 MHz, 
Figs. El E] and 00 show that our scheme requires measure- 
ment of a 1 pA current that changes in the nanosecond 
time scale. This is at the edge of the current technol- 
ogy that allows the measurement of 1 pA current with 
dynamics in the GHz frequency range with repeated- 
measurement techniqueiiiLiSii2iS£ 

One issue we have been trying to address in this study 
is to compare the measurement capability of a QPC de- 
tector and an SET detector. In terms of the theoretical 
descriptions of the qubit-detector interaction, the major 
difference between the QPC detector studied here and 
the SET detector in Refi& is that we model each QPC 
by a tunnel junction (Ref, 9 ), so that the QPC-qubit in- 
teraction directly modifies strength of tunneling, while in 
Ref£, the SET-qubit interaction influences both the SET 
island state energy and the island-lead tunneling. De- 
spite these differences, our numerical results showed that 
the current through the coupled QPC exhibits behaviors 
similar to those of the two-island SET current in many 
respects, such as in distinguishing the different qubit 
product states, in distinguishing the Bell-type entangled 
states from the product states, and in the observation 
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FIG. 6: Time-dependence of p AA (t),p BB (i), p cc ' (t) and 

p DD (t) for the 17 = case (rL ±5 ' = ri ±) ), starting from 
p DD (t = 0) = 1, for different measurement strengths (in 
terms of Td)- Here the intra-qubit tunneling rates are Q,l = 
Q,r — 0.75F, and there is no interaction between the qubits: 
J = 0. As measurement strength Td increases, the coherent 
motions of qubits slow down, which is a clear evidence of the 
quantum Zeno effect. 
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FIG. 7: Time dependent concurrence of a two-qubit state 
starting from a singlet state |e4) as a function of the dephasing 
rate Fd in the same parameter region as Fig llVI At t = 
and Td = 0, the concurrence takes a value of 1 and rapidly 
decreases to zero for large dephasing rates. 



of quantum Zeno effect for the qubit product states. 
Stronger differences between QPC and SET detectors do 
appear when the qubit-detector interaction strength in- 
creases. The measurement current of the detector that 
has a larger number of internal degrees of freedom (the 
two-island SET) seems to be able to describe more elab- 
orate quantum states of the two qubits. For example, 
the SET current can clearly distinguish the four product 
states shown in Ref 8 ; , while with the present QPC de- 
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FIG. 8: Time dependent currents for |e4) (singlet state) (panel 
(a)) and |ea) state (panel (b)), and the diagonal density ma- 
trix element (panel (c)) for the singlet state, when the de- 
phasing rate Td is changed with A = 0. The parameters are 
the same as those in Fig. At T d < 0.04F, Fig. (a) presents 
the proof of time- independence of the singlet state in Eg, 1191 1. 
and Fig.(b) shows the proof of the 4Q oscillation of Eg. 1201 . 
Figure (c) shows that the two-qubit state begins to include 
states other than the singlet states, resulting in the oscillation 
of the current (panel (a)) when Td becomes large. 

tector the current shows a simpler structure and smaller 
differences for the different qubit states. Qualitatively, 
the tunneling rate of a QPC is generally larger than that 
of an SET, which corresponds to shorter dwelling time 
for the QPC (in the present study the dwelling time for 
QPCs is effectively taken to be zero). This difference es- 
sentially originates from the simpler structure of a QPC 
compared to an SET. 

In the present study we obtained the density-matrix 
equations under the condition that the voltage bias be- 
tween the left and right electrodes is sufficiently large 
such that the left-right symmetry is broken and the trans- 
mission of electrons from the right electrode to the left 
can be neglected. Thus we cannot directly calculate the 
QPC differential conductance, which would provide more 



information for some experiments^ 1 - This is one of the 
limitations of the present method. An approach that can 
properly deal with low bias situations is still in develop- 
ment. 

Our configuration of the qubit-QPC coupling scheme 
can be straightforwardly extended to N(N > 2)-qubit 
detection. However, it depends strongly on the sensitiv- 
ity of the current readout circuit such that the 2 N states 
can be differentiated^ and is thus better suited for only 
a small number of qubits. In any case, the key objective 
of the present study is to obtain two-qubit information 
directly and dynamically, not to invent a general detector 
for a multi-qubit system, for which other configurations 
such as a typical one-detector-per-qubit setup are proba- 
bly more suitable and have to be further studied both ex- 
perimentally and theoreticallyii2ii2i22i2ii2£ Furthermore, 
we have considered an ideal measurement process in the 
present study. In a more realistic situation, imperfec- 
tions such as gate operation errors^ charge fluctuations 
around the qubit-QPC systems^ and phonons have to 
be considered. These imperfections could seriously re- 
duce the sensitivity of a measuring device. Thus more 
detailed analysis for the coupled multiqubit-detector sys- 
tem needs to be carried out in the future to further clarify 
these issues. 



VI. CONCLUSION 

We have solved master equations and described vari- 
ous time-dependent measurement processes of two charge 
qubits by two serially-coupled QPCs. The current 
through the QPCs is shown to be an effective means to 
measure the results of quantum calculations and the pres- 
ence of entangled states. Two-qubit dynamics is studied 
analytically and it is found that period of coherent os- 
cillation depends on their initial state. Our results thus 
show that the serially-coupled QPC can be an effective 
detector of two-qubit states of a pair of (coupled) charge 
qubits. 
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APPENDIX A: DERIVATION OF DENSITY MATRIX EQUATIONS 



Here we display all the density matrix equations of the qubit-QPC system. The density matrix equations can be 
classified according to the electronic states of the QPC island (See Fig. and qubit states (zi,Z2 — A,B,C.D, 

8=1,1). 



dt 



9 



dpi 1 ; 2 

dt 



dp* 12 
dt 



p(zi)-p(22)/ Z1Z2 | n ZlZ2\ 

1 R 1 R \Pb-\ "t" /\l / ) 



(Zl)'_ L p(»)'_ LT -l(zi)_ LF (^2) 



') 



(i[j X2 - J Z1 ]- [r<f y + r^ 2) ']) P ^ Z2 -iO fl (pf {Z1) ' Z2 -pl^^)-i^ L {p^^- P l^^) 



(Al) 



(A2) 



(A3) 



where 



r (A) 

1 L 
1 ft 
1 L 



r 



(A)' 



and 



_ r (B) 
L 

_ r (C) 

_ r (B)' 
L 



= r 
= r 



r 



(C)' 



C+) 

L ' 
(+) 
R ' 
r (+)' 

r (+)' 
1 R ' 



p (C) _ r (D) _ p(-) 
L — L — L > 
p(B) _ p (D) _ r (-) 



i?. 
r (C)' 
1 L 
r (B)' 
1 ft 



ft. 
r (D)' 

r (£>)' 

1 R 



ft 

r 



(-)' 

r (-)' 



gi(zi) and g r {zi) are introduced for the sake of notational 
convenience and represent relationships between different 
two-qubit states in the equations for the density matrix 
elements: 



g r (A) =B, 9l {A)=C, g r {B)=A, 9l (B) : 
g r (C) = D, 9l (C) = A, g r (D) = C :9l (D) 



D, 
-- B. 



J A =A L + A fl +J, J B =A L -A R -J, 
J c =-A L + A R -J, J D =~A L -A R +J. 
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